


	******************************************************************
	**
	**
	**		NAME:	ELITE PERSISTENCE IN THE ERA OF EXPANDING OVERSEAS TRADE
	**
	**		MAKES:  Appendix E Figures (a) and (b)		       
	**		
	**
	******************************************************************
	
	
	
	
	

			

		
			**	**	**	**	**	**	**	**	**	**	**	**	**
			**	**	**	**	**	**	**	**	**	**	**	**	**
			**	**	**	**	**	**	**	**	**	**	**	**	**
	


*-------------------------------------------------------------------------------
* global for input/output
*-------------------------------------------------------------------------------


* input
*------

global input = "" 	
	
	
* output
*-------

global output = ""	
		
	
	
			**	**	**	**	**	**	**	**	**	**	**	**	**
			**	**	**	**	**	**	**	**	**	**	**	**	**
			**	**	**	**	**	**	**	**	**	**	**	**	**
	


*-------------------------------------------------------------------------------
* preliminaries
*-------------------------------------------------------------------------------


* preliminaries
*--------------

clear
set more off





			**	**	**	**	**	**	**	**	**	**	**	**	**
			**	**	**	**	**	**	**	**	**	**	**	**	**
			**	**	**	**	**	**	**	**	**	**	**	**	**
	


*-------------------------------------------------------------------------------
* dataset
*-------------------------------------------------------------------------------



* use
*----

cd ${input}
use "analysis_vars_parl.dta", clear



			**	**	**	**	**	**	**	**	**	**	**	**	**
			**	**	**	**	**	**	**	**	**	**	**	**	**
			**	**	**	**	**	**	**	**	**	**	**	**	**
	


*-------------------------------------------------------------------------------
* assert
*-------------------------------------------------------------------------------


* assert
*-------

assert 	parl_year >= 1536 & parl_year <= 1760


duplicates tag parl_group, gen(tt)
assert tt == 0 if parl_group != .
drop tt






			**	**	**	**	**	**	**	**	**	**	**	**	**
			**	**	**	**	**	**	**	**	**	**	**	**	**
			**	**	**	**	**	**	**	**	**	**	**	**	**
	


*-------------------------------------------------------------------------------
* globals
*-------------------------------------------------------------------------------


* globals
*--------

global yr1 		= "parl_year > 1553 & parl_year < 1639"
global yr2 		= "parl_year >= 1660 & parl_year <= 1756"

global yr3 		= "parl_year >= 1639 & parl_year <= 1649"
global yr4 		= "parl_year >= 1628 & parl_year <= 1660"

global yr5 		= "parl_year < 1639 & parl_year >= 1583"
global yr6 		= "parl_year > 1550 & parl_year <= 1750"

global yr7 		= "parl_year > 1553 & parl_year <= 1640"


global xline1	= "xline(1660, lcolor(gs12))"
global xline2	= "xline(1688, lcolor(gs12) lpattern(line) )"
global xline3	= "xline(1640, lcolor(gs12))"
global xline4	= "xline(1562, lcolor(gs12)  lpattern(line) )"


	
	


			**	**	**	**	**	**	**	**	**	**	**	**	**
			**	**	**	**	**	**	**	**	**	**	**	**	**
			**	**	**	**	**	**	**	**	**	**	**	**	**
	

	
		
			**	**	**	**	**	**	**	**	**	**	**	**	**
			**	**	**	**	**	**	**	**	**	**	**	**	**
			**	**	**	**	**	**	**	**	**	**	**	**	**


*-------------------------------------------------------------------------------
* PLOT:  other commercial (FIGURE E1)
*-------------------------------------------------------------------------------




* hi lo
*------

capture drop hi 
capture drop lo
gen hi = 31
gen lo = 0
sort parl_group



* this plot
*----------

global outcome1 = "int_commtot_4_pperc"  

replace int_commtot_4_pperc = int_commtot_4_pperc * 100

assert int_commtot_4_pperc <= 100 & int_commtot_4_pperc >= 0 if int_commtot_4_pperc != .
global label1 = "% MPs Occupied in Commerce that is not Overseas Atlantic Trade"

  
global scale = "0(10)30"

  

* plot
*-----

# delimit ;

twoway 
	
	
		// shaded area
	(rarea lo hi parl_year	
		if   parl_year >= 1639 & parl_year <= 1649 
		,  color(gs12) )
		
		

	(lpoly ${outcome1} parl_year
		if   ${yr1}		
		,  bwidth(8) lcolor(gs1)  lwidth(medthick) lpattern(line)  ${xline2} ${xline4} )
	(lpoly ${outcome1} parl_year
		if   ${yr2}		
		,  bwidth(8) lcolor(gs1)  lwidth(medthick) lpattern(line)  )

		
		
		
	(scatter ${outcome1} parl_year
		if    ${yr1}		
		,  mcolor(gs1) msize(medlarge) msymbol(circle))
	(scatter ${outcome1} parl_year
		if    ${yr2}			
		,  mcolor(gs1) msize(medlarge) msymbol(circle))		


		
		
		,
		yscale(noline)
		xscale(noline)

		
		xsize(8.5)
		ysize(3)

		title("", color(black) size(medsmall) pos(11) ) 
		xlab(  1550(25)1750, 
			labsize(medsmall) nogrid) 
		ylab( ${scale}, 
			angle(hori) labsize(medsmall) nogrid  )
		xtitle("Year" , 
			color(black) size(large) )
		ytitle("Percentage of MPs", 
			color(black) size(large)  )
		graphregion( fcolor(white) lcolor(white) ) 
		plotregion( fcolor(white) lstyle(none) lcolor(white) ilstyle(none)) 
		legend(order(				
				2 	"${label1}"

				)
			cols(1)
			pos(6)
			
			region( color(none) )
			size(large)
			)
		
		;

# delimit cr


	* export graph
	*-------------

	cd "${output}"
	graph export "trend_commercial_othercomm.pdf", replace
	
	
	
	
			**	**	**	**	**	**	**	**	**	**	**	**	**
			**	**	**	**	**	**	**	**	**	**	**	**	**
			**	**	**	**	**	**	**	**	**	**	**	**	**


*-------------------------------------------------------------------------------
* PLOT:  trend in non commercial capital (FIGURE E2)
*-------------------------------------------------------------------------------


* hi lo
*------

capture drop hi lo
gen hi = 31
gen lo = 0
sort parl_group



* this plot
*----------

global outcome2 = "int_noncommcapprof_pperc"

replace int_noncommcapprof_pperc = int_noncommcapprof_pperc * 100 
assert int_noncommcapprof_pperc >= 0 & int_noncommcapprof_pperc <= 100 if int_noncommcapprof_pperc != .

global label = "all constituencies"
global scale = "0(10)30"



* plot
*-----

# delimit ;

twoway 
	
	
		// shaded area
	(rarea hi lo parl_year	
		if   ${yr3} 
		,  color(gs12) )


		
	(lpoly ${outcome2} parl_year
		if   ${yr1}		
		,  bwidth(7) lcolor(gs1)  lwidth(medthick) lpattern(line)  ${xline2} ${xline4} )
	(lpoly ${outcome2} parl_year
		if   ${yr2}		
		,  bwidth(7) lcolor(gs1)  lwidth(medthick) lpattern(line)  )
		
		
	(scatter ${outcome2} parl_year
		if    ${yr1}		
		,  mcolor(gs1) msize(medlarge) msymbol(circle))
	(scatter ${outcome2} parl_year
		if    ${yr2}			
		,  mcolor(gs1) msize(medlarge) msymbol(circle))		
		
		,
		yscale(noline)
		xscale(noline)
		
		
		xsize(8.5)
		ysize(3)

		title("", color(black) size(medsmall) pos(11) ) 
		xlab( 1550(25)1750, 
			labsize(medsmall) nogrid) 
		ylab( ${scale}, 
			angle(hori) labsize(medsmall) nogrid  )
		xtitle("Year" , 
			color(black) size(large) )
		ytitle("Percentage of MPs", 
			color(black) size(large)  )
		graphregion( fcolor(white) lcolor(white) ) 
		plotregion( fcolor(white) lstyle(none) lcolor(white) ilstyle(none)) 
		legend(order(
				2 "% MPs w/ Non-Commercial Capitalist Occupations"
				)
			cols(2)
			pos(6)
			
			region( color(none) )
			size(large)
			)
		
		;

# delimit cr


	* export graph
	*-------------

	cd "${output}"
	graph export "trend_capital.pdf", replace
	
	
	
			**	**	**	**	**	**	**	**	**	**	**	**	**
			**	**	**	**	**	**	**	**	**	**	**	**	**
			**	**	**	**	**	**	**	**	**	**	**	**	**
	


*-------------------------------------------------------------------------------
* PLOT:  trend in % pre-1550 family (FIGURE E3 (a))
*-------------------------------------------------------------------------------


* hi lo
*------

capture drop hi lo
gen hi = 101
gen lo = 0
sort parl_group



* this plot
*----------

global outcome1 = "ln_group_pre1550_n_eng_pperc"
replace ln_group_pre1550_n_eng_pperc = ln_group_pre1550_n_eng_pperc * 100
assert ln_group_pre1550_n_eng_pperc >= 0 & ln_group_pre1550_n_eng_pperc <= 100 if ln_group_pre1550_n_eng_pperc != .

global scale = "0(20)100"



* plot
*-----

# delimit ;

twoway 
	
	
		// shaded area
	(rarea hi lo parl_year	
		if   ${yr3} 
		,  color(gs12) )
		

		
	(lpoly ${outcome1} parl_year
		if   ${yr1}		
		,  bwidth(8) lcolor(gs1)  lwidth(medthick) lpattern(line)   )
	(lpoly ${outcome1} parl_year
		if   ${yr2}		
		,  bwidth(8) lcolor(gs1)  lwidth(medthick) lpattern(line) ${xline2} ${xline4} )	
		

		
		

		
	(scatter ${outcome1} parl_year
		if    ${yr1}		
		,  mcolor(gs1) msize(medlarge) msymbol(circle))
	(scatter ${outcome1} parl_year
		if    ${yr2}			
		,  mcolor(gs1) msize(medlarge) msymbol(circle))
		

		
		,
		yscale(noline)
		xscale(noline)

		
		xsize(8.5)
		ysize(3)

		title("", color(black) size(medsmall) pos(11) ) 
		xlab( , 
			labsize(medsmall) nogrid) 
		ylab( ${scale}, 
			angle(hori) labsize(medsmall) nogrid  )
		xtitle("Year" , 
			color(black) size(large) )
		ytitle("Percentage of MPs", 
			color(black) size(large)  )
		graphregion( fcolor(white) lcolor(white) ) 
		plotregion( fcolor(white) lstyle(none) lcolor(white) ilstyle(none)) 
		legend(order(
				2 "% MPs w/ a Surname From an MP in Parliament Before 1550"
				)
			cols(3)
			pos(6)
			
			region( color(none) )
			size(large)
			)
		
		;

# delimit cr


	* export graph
	*-------------

	cd "${output}"
	graph export "trend_pre1550.pdf", replace	
	
	
	
	
	
	
	
	
			**	**	**	**	**	**	**	**	**	**	**	**	**
			**	**	**	**	**	**	**	**	**	**	**	**	**
			**	**	**	**	**	**	**	**	**	**	**	**	**
	


*-------------------------------------------------------------------------------
* PLOT:  trend in commercial family (FIGURE E3 (b))
*-------------------------------------------------------------------------------


* hi lo
*------

capture drop hi lo
gen hi = 20
gen lo = 0
sort parl_group



* this plot
*----------

global outcome1 = "fam_any_comm_pperc"
replace fam_any_comm_pperc = fam_any_comm_pperc * 100
assert fam_any_comm_pperc >= 0 & fam_any_comm_pperc <= 100 if fam_any_comm_pperc != .

global scale = "0(5)20"



* plot
*-----

# delimit ;

twoway 
	
	
		// shaded area
	(rarea hi lo parl_year	
		if   ${yr3} 
		,  color(gs12) )
		

		
	(lpoly ${outcome1} parl_year
		if   ${yr1}		
		,  bwidth(8) lcolor(gs1)  lwidth(medthick) lpattern(line)   )
	(lpoly ${outcome1} parl_year
		if   ${yr2}		
		,  bwidth(8) lcolor(gs1)  lwidth(medthick) lpattern(line) ${xline2} ${xline4} )	
		

		
		

		
	(scatter ${outcome1} parl_year
		if    ${yr1}		
		,  mcolor(gs1) msize(medlarge) msymbol(cicle))
	(scatter ${outcome1} parl_year
		if    ${yr2}			
		,  mcolor(gs1) msize(medlarge) msymbol(cicle))
		

		
		,
		yscale(noline)
		xscale(noline)
		
		
		xsize(8.5)
		ysize(3)

		title("", color(black) size(medsmall) pos(11) ) 
		xlab( , 
			labsize(medsmall) nogrid) 
		ylab( ${scale}, 
			angle(hori) labsize(medsmall) nogrid  )
		xtitle("Year" , 
			color(black) size(large) )
		ytitle("Percentage of MPs", 
			color(black) size(large)  )
		graphregion( fcolor(white) lcolor(white) ) 
		plotregion( fcolor(white) lstyle(none) lcolor(white) ilstyle(none)) 
		legend(order(
				2 "% MPs w/ Family Occupied in Commerce"
				)
			cols(3)
			pos(6)
			
			region( color(none) )
			size(large)
			)
		
		;

# delimit cr


	* export graph
	*-------------

	cd "${output}"
	graph export "trend_commfamily.pdf", replace
	
	
	
	
	
			**	**	**	**	**	**	**	**	**	**	**	**	**
			**	**	**	**	**	**	**	**	**	**	**	**	**
			**	**	**	**	**	**	**	**	**	**	**	**	**
	


*-------------------------------------------------------------------------------
* PLOT: average dynasty age (FIGURE E4 (a))
*-------------------------------------------------------------------------------
	
	

* hi lo
*------

capture drop hi lo
gen hi = 300
gen lo = 100
sort parl_group



* this plot
*----------

global outcome1 = "dynas_age_nation_avg" 
global scale = "100(50)300"



* plot
*-----

# delimit ;

twoway 
	
	
		// shaded area
	(rarea hi lo parl_year	
		if   ${yr3} 
		,  color(gs12) )
		

		
	(lpoly ${outcome1} parl_year
		if   ${yr1}		
		,  bwidth(8) lcolor(gs1)  lwidth(medthick) lpattern(line)   )
	(lpoly ${outcome1} parl_year
		if   ${yr2}		
		,  bwidth(8) lcolor(gs1)  lwidth(medthick) lpattern(line) ${xline2} ${xline4} )	
			
		
		

		
	(scatter ${outcome1} parl_year
		if    ${yr1}		
		,  mcolor(gs1) msize(medlarge) msymbol(circle))
	(scatter ${outcome1} parl_year
		if    ${yr2}			
		,  mcolor(gs1) msize(medlarge) msymbol(circle))
		
		
		,
		yscale(noline)
		xscale(noline)
		

		
		xsize(8.5)
		ysize(3)

		title("", color(black) size(medsmall) pos(11) ) 
		xlab( , 
			labsize(medsmall) nogrid) 
		ylab( ${scale}, 
			angle(hori) labsize(medsmall) nogrid  )
		xtitle("Year" , 
			color(black) size(large) )
		ytitle("Years since first surname appearance", 
			color(black) size(large)  )
		graphregion( fcolor(white) lcolor(white) ) 
		plotregion( fcolor(white) lstyle(none) lcolor(white) ilstyle(none)) 
		legend(order(
				2 "Average Years Since First Surname Appearance for MPs w/ a Pre-Existing Surname"
				)
			cols(3)
			pos(6)
			
			region( color(none) )
			size(large)
			)
		
		;

# delimit cr


	* export graph
	*-------------

	cd "${output}"
	graph export "trend_avgdynasage.pdf", replace	
	
	
	
	

			**	**	**	**	**	**	**	**	**	**	**	**	**
			**	**	**	**	**	**	**	**	**	**	**	**	**
			**	**	**	**	**	**	**	**	**	**	**	**	**
	


*-------------------------------------------------------------------------------
* PLOT: new dynasties being started (FIGURE E4 (b))
*-------------------------------------------------------------------------------
	
	
* hi lo
*------

capture drop hi lo
gen hi = 100
gen lo = 0
sort parl_group



* this plot
*----------

global outcome1 = "dynas_first_nation_pperc"  
replace dynas_first_nation_pperc = dynas_first_nation_pperc * 100 
assert dynas_first_nation_pperc >= 0 & dynas_first_nation_pperc <= 100 if dynas_first_nation_pperc != . 
global scale = "0(20)100"



* plot
*-----

# delimit ;

twoway 
	
	
		// shaded area
	(rarea hi lo parl_year	
		if   ${yr3} 
		,  color(gs12) )
		

		
	(lpoly ${outcome1} parl_year
		if   ${yr1}		
		,  bwidth(8) lcolor(gs1)  lwidth(medthick) lpattern(line)   )
	(lpoly ${outcome1} parl_year
		if   ${yr2}		
		,  bwidth(8) lcolor(gs1)  lwidth(medthick) lpattern(line) ${xline2} ${xline4} )	
			
		
		

		
	(scatter ${outcome1} parl_year
		if    ${yr1}		
		,  mcolor(gs1) msize(medlarge) msymbol(circle))
	(scatter ${outcome1} parl_year
		if    ${yr2}			
		,  mcolor(gs1) msize(medlarge) msymbol(circle))
		
		
		,
		yscale(noline)
		xscale(noline)
		
		
		xsize(8.5)
		ysize(3)

		title("", color(black) size(medsmall) pos(11) ) 
		xlab( , 
			labsize(medsmall) nogrid) 
		ylab( ${scale}, 
			angle(hori) labsize(medsmall) nogrid  )
		xtitle("Year" , 
			color(black) size(large) )
		ytitle(" ", 
			color(black) size(large)  )
		graphregion( fcolor(white) lcolor(white) ) 
		plotregion( fcolor(white) lstyle(none) lcolor(white) ilstyle(none)) 
		legend(order(
				2 "% MPs w/ a New-to-Politics Surnames in the First Parliament of their Appearance"
				)
			cols(3)
			pos(6)
			
			region( color(none) )
			size(large)
			)
		
		;

# delimit cr


	* export graph
	*-------------

	cd "${output}"
	graph export "trend_percnewdynas.pdf", replace	
	
	
		

			**	**	**	**	**	**	**	**	**	**	**	**	**
			**	**	**	**	**	**	**	**	**	**	**	**	**
			**	**	**	**	**	**	**	**	**	**	**	**	**
	


*-------------------------------------------------------------------------------
* PLOT:  trend in dynastic (england only) (FIGURE E5 (a))
*-------------------------------------------------------------------------------


* hi lo
*------

capture drop hi 
capture drop lo
gen hi = 10
gen lo = .4
sort parl_group



* this plot
*----------

global outcome1 = "dynas_all_nationE_pperc" 
replace dynas_all_nationE_pperc = dynas_all_nationE_pperc * 100 
assert dynas_all_nationE_pperc >= 0 & dynas_all_nationE_pperc <= 100 if dynas_all_nationE_pperc != .
global scale = "0(20)100"



* plot
*-----

# delimit ;

twoway 
	
	
		// shaded area
	(rarea hi lo parl_year	
		if   ${yr3} 
		,  color(gs12) )
		
		

	(lpoly ${outcome1} parl_year
		if   ${yr1}		
		,  bwidth(8) lcolor(gs1)  lwidth(medthick) lpattern(line)  ${xline2} ${xline4} )
	(lpoly ${outcome1} parl_year
		if   ${yr2}		
		,  bwidth(8) lcolor(gs1)  lwidth(medthick) lpattern(line)  )
		
		
		
		
		
	(scatter ${outcome1} parl_year
		if    ${yr1}		
		,  mcolor(gs1) msize(medlarge) msymbol(circle))
	(scatter ${outcome1} parl_year
		if    ${yr2}			
		,  mcolor(gs1) msize(medlarge) msymbol(circle))		

		
		,
		yscale(noline)
		xscale(noline)
		

		
		xsize(8.5)
		ysize(3)

		title("", color(black) size(medsmall) pos(11) ) 
		xlab(  1550(25)1750, 
			labsize(medsmall) nogrid) 
		ylab( ${scale}, 
			angle(hori) labsize(medsmall) nogrid  )
		xtitle("Year" , 
			color(black) size(large) )
		ytitle("Percentage of MPs", 
			color(black) size(large)  )
		graphregion( fcolor(white) lcolor(white) ) 
		plotregion( fcolor(white) lstyle(none) lcolor(white) ilstyle(none)) 
		legend(order(
				2 "% MPs w/ Pre-existing Surnames (England-only)"
				)
			cols(3)
			pos(6)
			
			region( color(none) )
			size(large)
			)
		
		;

# delimit cr


	* export graph
	*-------------

	cd "${output}"
	graph export "trend_dynas_england.pdf", replace	


	
			**	**	**	**	**	**	**	**	**	**	**	**	**
			**	**	**	**	**	**	**	**	**	**	**	**	**
			**	**	**	**	**	**	**	**	**	**	**	**	**
	


*-------------------------------------------------------------------------------
* PLOT:  trend in dynastic (wales only) (FIGURE E5 (b))
*-------------------------------------------------------------------------------


* hi lo
*------

capture drop hi 
capture drop lo
gen hi = 101
gen lo = 0
sort parl_group



* this plot
*----------

global outcome1 = "dynas_all_nationW_pperc" 
replace dynas_all_nationW_pperc = dynas_all_nationW_pperc * 100 
assert dynas_all_nationW_pperc >= 0 & dynas_all_nationW_pperc <= 100 if dynas_all_nationW_pperc != .
global scale = "0(20)100"



* plot
*-----

# delimit ;

twoway 
	
	
		// shaded area
	(rarea hi lo parl_year	
		if   ${yr3} 
		,  color(gs12) )
		
		

	(lpoly ${outcome1} parl_year
		if   ${yr1}		
		,  bwidth(8) lcolor(gs1)  lwidth(medthick) lpattern(line)  ${xline2} ${xline4} )
	(lpoly ${outcome1} parl_year
		if   ${yr2}		
		,  bwidth(8) lcolor(gs1)  lwidth(medthick) lpattern(line)  )
		
		
		
		
		
	(scatter ${outcome1} parl_year
		if    ${yr1}		
		,  mcolor(gs1) msize(medlarge) msymbol(circle))
	(scatter ${outcome1} parl_year
		if    ${yr2}			
		,  mcolor(gs1) msize(medlarge) msymbol(circle))		

		
		,
		yscale(noline)
		xscale(noline)
		

		
		xsize(8.5)
		ysize(3)

		title("", color(black) size(medsmall) pos(11) ) 
		xlab(  1550(25)1750, 
			labsize(medsmall) nogrid) 
		ylab( ${scale}, 
			angle(hori) labsize(medsmall) nogrid  )
		xtitle("Year" , 
			color(black) size(large) )
		ytitle("Percentage of MPs", 
			color(black) size(large)  )
		graphregion( fcolor(white) lcolor(white) ) 
		plotregion( fcolor(white) lstyle(none) lcolor(white) ilstyle(none)) 
		legend(order(
				2 "% MPs w/ Pre-existing Surnames (Wales-only)"
				)
			cols(3)
			pos(6)
			
			region( color(none) )
			size(large)
			)
		
		;

# delimit cr


	* export graph
	*-------------

	cd "${output}"
	graph export "trend_dynas_wales.pdf", replace	
	
	

	
	
	




			**	**	**	**	**	**	**	**	**	**	**	**	**
			**	**	**	**	**	**	**	**	**	**	**	**	**
			**	**	**	**	**	**	**	**	**	**	**	**	**
	

						** end of do file ** 	
				